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Abstract 

(N ■ 

^ I We present some dynamic and entropic considerations about the evolution of a continuous 

t^^ ' time quantum walk implementing the clock of an autonomous machine. On a simple model, 

(<— ^ . we study in quite explicit terms the Lindblad evolution of the clocked subsystem, relating 

tJ- I the evolution of its entropy to the spreading of the wave packet of the clock. We explore 

^ ■ possible ways of reducing the generation of entropy in the clocked subsystem, as it amounts 

Q . to a deficit in the probability of finding the target state of the computation. We are thus 

'^ I lead to examine the benefits of abandoning some classical prejudice about how a clocking 

pi i' mechanism should operate. 
-(— > 

c, 

§ ■ 1 Introduction 



>< 



We study a model of quantum computation originally due to Feynman ^lOj, in which the evolution 
of the quantum register is controlled by an auxiliary clock or cursor register. Feynman pointed 



;h ' out the computational relevance, for such a system, of a Hamiltonian of the form 

H = --J2 U{x,y;(T)T+iy)T4x) + h.c. (1) 

x,yeG 

by showing that the CCNOT primitive (and therefore a universal quantum computer) can be 
implemented by a suitable choice of the graph G and of the dependence on the register spins cr of 
the dynamical variables U which couple the cursor spin creation and annihilation operators t±. 
A general architecture emerges from the above model in which an excitation of the r field performs 
a quantum walk on a planar graph, the basic events being the flipping of a a spin determined by 
the flowing of the cursor current (the NOT primitive) and the conditional choice, determined by 
the state of a a qubit, of the flow of the cursor along alternative edges of the graph (the SWITCH 
primitive) . 

It is to the clock subsystem that we devote most of our attention in this paper, contributing an 
explicitly solvable toy model to the long-standing exploration j2Zj of the quantum limitations of 
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time measurement (for a discussion of the foundational nature of this problem and for extensive 
updated references, see |I2])- We study, furthermore, by frankly heuristic numerical means when 
necessary, variations on the theme of our toy model, with the purpose of gaining and sharing 
experience on a problem (optimizing the performance of a quantum clock) that might become of 
practical relevance in studying ballistic computation on nanostructures (3tlT7j. 
The dynamics under the Hamiltonian (^ was examined by A. Peres [SE! on a linear chain with 
edges between nearest neighbor sites, with the f/'s taken to be numerical functions of the positions; 
equivalently stated, the dynamics was studied along the linear chain of logical successors of an 
initial register state. Peres showed, in the suggestive terminology of the more recent Reference j7] , 
that unit fidelity could be achieved in the transfer of a state along such a linear chain by suitable 
engineering of the coupling constants U. 

Interest in the case in which the f/'s are numerical functions (as opposed to matrices acting on 
additional spins) has been revived by the success of the quantum walk paradigm j5j. In |7j the 
graph is considered in its all important role as a wire able to spatially transfer a single-spin quantum 
state; in ^ the interesting possibility is examined of performing the mirror inversion of a many- 
spin state. 

For a recent presentation of the full model (^, which we call an interacting XY system, we refer 
to jH] , where particular attention is paid to the role of additional controlling spins in implementing 
successive visits to selected parts of the graph in iterated computations (quantum subroutines) 
and in the storage of results in telomeric chains. 

In this note we examine an explicitly solvable, yet computationally non trivial, instance of the 
Hamiltonian (^, paying specific attention to non-positional observables of the system: our main 
concerns will be speed (of computation) and entropy (of controlled and/or controlling subsystem). 
We also consider the observable number of particles (agents performing a quantum walk along the 
XY chain), and discuss the interest and limitations of the proposal of a multi-hand quantum clock 
as a substitute for the loops implementing iterated applications of quantum subroutines. 

The paper is organized as follows. In Section 2 we present the model and establish our notation. 
In Section 3 we relate the speed of computation to the group velocity of the motion of the cursor 
wave packet along the graph. In Section 4 we discuss the build-up of entropy in the clocked 
subsystem because of the spreading of the wave packet of the clocking agent. An outline of 
possible choices of the initial form of this wave packet bringing such entropy build-up close to 
a minimum is given in Section 5. Section 6 examines the possible interest of multiagent spin 
networks. Section 7 is devoted to conclusions and outlook. 

2 The model 

The model we consider consists of two functionally distinct parts, the input/output register and 
the clock or cursor, and evolves as an autonomous system under the sole effect of its initial state 
not being an eigenstate of the Hamiltonian of the composite system. 
Let I i?(l) ) G Hregister be the initial state of the input/output register. Set 

d = dim(Hregister)- (2) 

Let Ui, U2, . . . , f^AT-i be the unitary operators representing the successive primitive steps of the 
computation to be performed. Suppose, namely, that the goal of the computation is to transform 
the input state | -R(l) ) into the output state | R{N) ) = U^-i ■ . . .U2 ■ Ui\ -R(l) ) by visiting the 



successive intermediate states 

I R(x) ) = f/^_i| i?(x - 1) ), 1 < a; < N. (3) 

Following the approach of PO], we model the clocking mechanism, which sequentially applies the 
transformations Ui,U2, ■ ■ ■ , Un^i to the register, with a quantum mechanical system, the cursor. 
We call Ticursor the s-dimensional {s > N) state space of this system and refer it to a selected 
orthonormal basis | C(l) ), | C(2) ),..., | C{s) ). It will help the intuition, and will explain the 
notation used below, to think of an explicit implementation of the cursor by s spin 1/2 particles 
and to think of the state \ C{x)) as obtained by flipping "up" the spin in position x with respect 
to the "all down" reference state. 

We suppose that the state of the overall system, the machine, evolves in the Hilbert space 
'Hjnachine = 'Hregister ® 'Hcursor Under the actiou of a Hamiltoniau of the form 

A "~^ 

F = -- J] f/. ® I C(x + 1) ){C{x) I + U-^ ® I C{x) ){C{x + l)\. (4) 



x=l 



Notice that only [/i, U2, ■ ■ ■ , Un-i are assigned by the algorithm we are interested in; Un, • • • , Us-i 

are to be assigned as a part of the description of the clocking mechanism. 

For instance. Reference |2j presents the case in which Un, • • • , Us-i = Ir, the identity in Hregister , 

and shows the role of the cursor sites A^, . . . , s as a storage mechanism of the output | R{N) ). 

An alternative point of view was taken in some of the numerical examples of |H^, motivated by 

Grover's algorithm: one may suppose all of the U^ to coincide, in such a context, with Grover's 

estimation ■ oracle step G, and study the effect of applying G more than the optimal number A^ — 1 

of times. 

This point of view will be taken also in some numerical examples of this paper, where we focus 

our attention on a Hamiltonian of the form: 

A "'^ 

H = --^f/^®r+(x + l)r_(x) + f/"^ ® r+(x)r_(x + 1). (5) 

x=l 

Most of our numerical examples will refer in fact to the following particular instance (Toy model): 

A '"^ 
Ht = --^e-*t'^2^+(x + l)r_(x) + e'?'"2r+(a;)r_(a; + l). (6) 

^ x=l 

Here a= (cti, cr2, o's) is the single register spin 1/2 that we are going to consider in our model; the 
cursor subsystem is implemented as a collection of spin 1/2 systems r(j) = {ti{j) , T2{j) , r^^j)) , 
j e 1, . . . , s, and T±{j) = (ri(j) ± iT2{j))/2. 

The Hamiltonian (JH) can, of course, be considered as the restriction of (jSj) to the eigenspace 
belonging to the eigenvalue 1 of 

N, = ±'-±^, (7) 

x=l 

provided we identify | C{x) ) with the simultaneous eigenstate of r3(l), r3(2) . . . , r3(s) in which only 
Tsi^x) has eigenvalue +1. 



3 Speed of computation 

The eigenvalue problem for the Haniiltonian (j3)) is solved by the following Ansatz for the eigenstates 



s-l 



I e ) = J2v{x) ?7,_i • . . . ■ f/i| r ) ® I C{x) ), (8) 

x=l 

with I r ) any non vanishing vector in Tiregister , suggested by the conservation laws discussed in 
j26j (equation 19, p. 3270). Inserting this Ansatz, the eigenvalue problem becomes: 

ev{x) = --{v{x+l)+v{x-l)). (9) 

It is immediate to recognize in the right hand side a finite difference approximation of the Laplace 
operator (the free Schrodinger equation is, not surprisingly, at work in the motion of the clock); 
together with the boundary conditions v{0) = v{s + 1) = 0, this leads in an obvious way to the 
eigenvalues 

Cfc = -A cos ( — — j , A; = 1, 2, . . . , s. (10) 

The multiplicity of each eigenvalue is equal to d = dim{l-iregister )• An orthonormal basis in the 
eigenspace belonging to the eigenvalue Cfc is given by: 

s-l 

I efc;r-j ) = '^Vk{x)U^-i ■ ...■Ui\rj) ®\ C{x) ), (11) 

x=l 

where | ri ),..., | r^ ) is an orthonormal basis in Hregister , and 



.,(x) = y^sin(-^xj. (12) 

The same statements hold, of course, for the Haniiltonian (0) in the eigenspace belonging to the 

eigenvalue 1 of A^s. 

An initial state (at time t = 0) of the form 

\M^) = \R{l))®\Cil)) (13) 

evolves, under the Haniiltonian (0]), into 

s 

I Ml (t) ) = e-'^* I Ml ) = ^ c(t, x; s) \ R{x) )®\C{x)) (14) 

x=l 

where 



s + l 



(15) 



c(t, x\ s) = > exp iXt cos ( I sin ( | sin ( 

s + 1^ [ \s + l)\ \s + l) V 

Equation (fT^ singles out the interest of the observable 

Q = Y.x I C{x) ){C{x) I = X^xi±^, (16) 

x=l x=l 



the position of the cursor: it shows that if, at any time t, Q is measured on a system in the state 
I Mi(t) ) and the the resuh x is observed, then the register collapses into the state | R{x) ) obtained 
from I -R(l) ) by the application of Ui, U2, ■ ■ ■ , U^-i in the right order. 

The observable Q/t acquires thus the meaning of number of primitives per unit time applied to 
the initial condition | -R(l) ) in the time interval (0,t). In order to study the behavior over long 
intervals of time {t -^ +cxd) of this observable in the case of a long computation (s -^ +00) it is 
expedient to study its characteristic function 



^^^ j I Mi(t) ) = ^ |c(t, x; s)\^exp (^izj^ 

^ ^ x=l 



0,,i(z) = (Mi(t) |exp (iz^]\ Mi{t) ) = V|c(t,x;s)|2exp (iz-) , (17) 



namely the Fourier transform of its probability distribution. 

The large s behavior is easily studied by inserting the explicit integral representation of the s — > 
+CXD limit of (fT^ into (jTTj) ; the t — » +00 limit is similarly studied by substituting the sum over 
V = x/t, step 1/t, appearing in (fTTj) with an integral and evaluating the leading contributions by 
a standard stationary phase argument. We thus obtain: 

lim lim ( Mi(t) I exp ( iz-] | Mi(t) ) = / — ^'" e^^dv. (18) 

t^oos-+oo \ tj Jo vrv 1 — f 2 

As convergence in the sense of characteric functions implies convergence in the sense of cumulative 
distribution functions {convergence in law), we conclude that a "long" computation starting from 
the initial condition ()1H|1 proceeds "in the long run" at a rate of V{Mi) steps per unit time (the 
unit of time having been set so that A = 1 ), V{Mi) being the random variable defined by having as 
its characteristic function the right hand side of (fTHj) : equivalently stated it has probability density 

/v(Mi)(f) = Ao,i)(^) — /- ^ (19) 

Here and in what follows we denote by /(a,6)the indicator function of an interval (a, b): 



The mean value 



E{V{M,)) = I V fviM,){v)dv = — (21) 



'0 

and the variance 

var{V{M,)) = E {{V{M,)y) - E {V{M,)f = ^ - (|^) ' (22) 

are then easy to compute from (fTn|) . 

More generally, for any positive integer xq, a state such as 

I M,„ ) = I R{x^) ) ® I C(a;o) ) (23) 

having at a certain instant the cursor in xq, evolves with a speed l^(Mj,.y) having cumulative 
distribution function 

FviM..M ^ ProbiViM,,)<v)= (24) 

, , /2arcsin(f) sin(2xoarcsin(t>))\ ^ , , 



and expectation value 

E(V(M,,)) = ^^^ (25) 

Comparison between (|2H) and (|25|) shows the effect of a measurement of Q. If, at a given t, Q 
is measured and the result xq is found, then the state (fT^ . into which the initial condition (fT^ 
has evolved, collapses into the state (J2S1)- From this moment on the computation proceeds at the 
mean rate ()25|) : for large values of t, reading the clock is likely to reduce the speed of further 
computation by a factor 3/4 (without, because of (fTT|) . altering its correctness). 

4 Entropy 

Motivated by the experience gained under the particular initial conditions (fT^ and (j^^ we define, 
for any (unentangled) initial condition of the form (for fixed e > 1) 

\R-i:o) = \R)^J2M^)\C{x)), (26) 

x=l 

the "time-of-fiight speed" |T1] of computation in the state ipo as the random variable V{iIjo) having 
characteristic function 

0v(^„)(z)= lim lim (i?;^o|e'*^expf2^^)e-^*^|i?;V'o). (27) 

The above limit is easily shown to exist by the techniques outlined in the previous section; it 
corresponds to the probability density 

|^(arcsin(t;))p + [^(tt - aTcsm{v))\^ 
/v(Vo)l^J = ^(0,1)1^) ;-, ., (28) 

where 

*b) = y~Yl sin(pa;)^o(a;). (29) 

The observable Q retains in this context the meaning of relational time ^2] in the sense that, given 

that at any parameter time t the cursor is found at x, it is then certain that the register is found 

in the state Ux-i ■ . . . ■ U2 ■ Ui\ R). 

In reading the output at any time t, namely in the measurement of any, however carefully chosen, 

observable of the register, there is an intrinsic uncertainty corresponding to the uncertainty about 

how far the computation has proceeded. The fact that Q/t has a non trivial limit in law means 

that the leading term of the variance of Q is proportional to t^ and therefore that the uncertainty 

increases with t. This section is devoted to the examination of an example in which the notion of 

"the most careful choice" of the observable to read on the register can be made precise and shown 

to be pertinent to the algorithm considered. 

We consider for the moment the initial condition | Mi ) given in ()13j) and its time evolution | Mi(t) ) 

described in (I14|l . More general initial conditions of the form ()2(')j) will be examined in the next 

section. 

Call 

p^(t) = |Mi(t))(Mi(t)| (30) 

6 



the density matrix of the machine at time t. 

By taking the partial trace Tr-j-c^^^.^^^ (pmit)) with respect to the cursor degrees of freedom, we get 

the density matrix Pr{t) of the register: 

s 

Pr{t) = Y, Ic(^'^5 s)\^\ R{x) ){R{x) \. (31) 

x=l 

Call \j(t) the positive eigenvalues of pr(t) and | bj(t) ) the corresponding eigenstates. A simple 
computation, amounting to the Schmidt decomposition [22] of the state (fT^ . shows, then, that 
the density matrix of the cursor is given by 

Pc(t) = X]A,(t)|rf,(t))(d,(t)| (32) 

j 

where 

I d,{t) ) = ^^ ^c(t,x; s){ h,{t) I R{x) ) I C{x) ). (33) 

Because of ()32p and of the orthonormality of the states | dj{t) ), the von Neumann entropy of the 
register and also of the cursor is then given by 

S{Pcit)) = -^A,(t)lnA,(t) = S{pr{t)). (34) 

i 

We observe that, as (|!?T|l shows, the von Neumann entropy of each subsystem does depend on the 

algorithm being performed. It is, indeed, only under the hypothesis, nowhere made above, that the 

states I R{x) ) are orthonormal that (jHTj) is the spectral decomposition of Pr{t) (the von Neumann 

entropy becoming in this case equal to the Shannon entropy of the distribution of Q). 

We focus our attention, in what follows, on our Toy model ©, in which the register is a single 

spin 1/2 system. We indicate by 61,62,63 the versors of the three coordinate axes to which the 

components g_= (ai, (72, era) of such a spin are referred. 

In the basis | 0-3 = ±1 ), the density operator Pr{t) will be represented by the matrix 

P'^^^- 2\s,{t)+is,{t) l-.3(t) ) ^^^^ 

where 

s,{t)=Tr{pr{t)-a,), J = 1,2,3. (36) 

Equivalently stated, the Block representative of the state pT{t) is given by the three-dimensional 
real vector 

s 

s_{t) = ^|c(t,x;s)|^(i?(x) \a\R{x)). (37) 



x=l 



We shall assume, in what follows, that the initial state of the cursor is | C(l) ) and that the initial 
state of the register is of the form 

I i?(l) ) = cos (^ I a3 = +1 ) + sin (^ ^3 = -1 ) (38) 



namely the eigenstate belonging to the eigenvalue +1 of n(l) ■ g_, with 

n{l) = CisinO + e^cos9. (39) 

We wish to remark that the above example captures the geometric aspects not only of such simple 
computational tasks as NOT or yNOT (viewed as rotations of an angle vr or 7r/2 respectively, 
decomposed into smaller steps of amplitude a) but also of Grover's quantum search |^. If, indeed, 
the positive integer yU is the length of the marked binary word to be retrieved, and we set 

x(/i) = arcsin(2"2 ) (40) 

and 

9 = 71-2 x(/i) (41) 

then the state (jHH|) correctly describes the initial state \ l) oi the quantum search as having a 
component 2~^/^ in the direction of the target state, here indicated by | cu ) = | as = +1 ), and a 
component \/l — 2^^ in the direction of the fiat superposition, here indicated by \ o^ = — 1 ), of 
the 2^ — 1 basis vectors orthogonal to the target state. In this notations, if 

a = -4x(/i), (42) 

then the unitary transformation exp(— zacr2/2) corresponds to the product B ■ A of the oracle step 

A = Ir-2\Lu){to\ (43) 

and the estimation step 

B = 2\l){l\- Ir. (44) 

We refer the reader to the beautifully pedagogical approach of Jozsa ^H] where it is shown that 

in Grover's search the //-qubits register evolves in the two dimensional space spanned by the its 

initial state | l ) and the target state \uj). Thus, one qubit suffices to represent all instances of 

quantum search. 

It is having in mind the connection with Grover's algorithm that, for the sake of definiteness, in the 

examples that follow we are going to consider the one-parameter family of models, parametrized 

by the positive integers fi, corresponding to the choice (j41|) and (j42p of the parameters 6 and a 

and to the choice s = 2^^ + 1 of the number of cursor sites, corresponding to the possibihty of 

performing up to an exhaustive search. 

In the example defined by the above conditions it is 

( R{x) I a\ R{x) ) = sin {6 + {x - l)a) e^ + cos {6 + {x - l)a) Cg (45) 

and, therefore, 

s 

s{t) = Y^ \c{t, x; s) f (sin {9 + {x - l)a) e^ + cos {9 + {x - l)a) Cg) . (46) 

x=l 

Figure 1 presents, inscribed in the unit circle, a parametric plot of {si{t) , s^{t)) under the above 
assumptions . It is convenient to describe the Bloch vector s_{t) = si(t) e^ + Ss(t) e^ in polar 
coordinates as 

si(t) = r(t) sin7(t), S3(t) = r{t) cos^{t). (47) 



% 




Sl 



Figure 1: A parametric plot of (si(t), S3(t)) for < t < s, A = 1. The choice /i = 7, x = 
arcsin(l/2'^''^), s = 2^^ + 1, a = —4%, ^ = vr — 2% of the parameters is motivated by the connection 
with Grover's algorithm. Only the initial state lies on the unit circumference, the locus of pure 

states. 

A very simple approximate representation of s_{t) becomes then possible: 



r{t)e 



ijit) 



y |c(t, x; s)p exp{i{6 + (x — l)a)) 



x=l 



exp{i{6 — a)) >^ |c(t, x; s)|^exp( 



zax) 



x=l 



exp{i{e - a))E {exp{iaXtV{Mi{0)))) 



(4J 



The last step, legitimate for 1 << A t < s, requires only the explicit computation of the charac- 
teristic function corresponding to the probability density (HD), which leads to 



r(t)e-«-^^^^P^^i^^^((J, 



(T) - T ,h{T) + i{TH,{T) - H,{T)))) 



(49) 



where Jk and H^ are, respectively, Bessel functions and Struve functions jSI], and T = aXt. 
The time evolution of the register subsystem is summarized by the Lindblad equation [T31 1201 



dt 



i d'lit) r , ,-, 1 (ilnr(t) , , , ,,, 

^^ ^ [a2,Pr{t)] + T— -^ [0^2, [a2,P.(t)]] . 



2 dt 



4 dt 



(50) 



The commutator term [cr2,Pr(^)] describes the Hamiltonian part of the dynamics (after all we are 
considering a rotation about the X2 axis); the double commutator [cr2, [cr2, Pr(^)]] describes, in much 
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S(Pr(t)) 

ln(2) 




150 200 250 

Figure 2: The von Neumann entropy of the register as a function of time, for the same model as 
in figure ^ for < t < s (sohd hne) and for s <t <2s (dashed hne). 



the same sense as equation 2.8 of [22], the decohering effect of this rotation being administered by 
the cursor in discrete steps at random times. 
The eigenvalues of pr{t) can be written as 

Ai(t) = l(l + r(t)), A2(t) = |(l-r(t)). (51) 

The von Neumann entropy S {pr(t)) is therefore 

5(p.W)^-i^l.i^-i^.ni^. (52) 

An example of its behaviour is shown in figure |21 The eigenvectors corresponding to the eigenvalues 
^^ are, respectively 

iMO)^(rS). iMO)^(;-<]<;)/^'). (53) 

It is to be stressed that, at each time t, the projector | bi(t) )(&i(i) | is, among the projectors on the 
state space of the register, the one having in the state pr(t) the greatest probability of assuming, 
under measurement, the value 1. Thus, the most careful choice (the one affected by minimum 
uncertainty) of the observable to read on the register at time t is the projector | bi{t) ){bi{t) \. In 
the case of Grover's algorithm one must measure the projector \uj){uj\ = \a3 = +1 )( cts = +1 | 
(and one easily can, because of the kickback mechanism analyzed, for instance, in jH]) and has the 
freedom of choosing the time r at which to perform the measurement. The best choice is therefore 
such that I 6i(r) ) = | a; ) (in our notational setting, r is the time at which the helix of figure ^ 
crosses for the first time the positive S3 axis). In spite of the fact of being now in the most favorable 
setting, one has, nevertheless, a deficit 1 — Ai(r) in the probability of finding the target state. 
As figure IHl shows, there are successive instants of time at which the probability of successful re- 
trieval has a local maximum (a remnant of the periodic nature of Grover's algorithm when applied 
by an outside macroscopic agent) but the heights of these successive maxima form a sequence 
having a decreasing trend. 
Further insight into our toy model is gained by examining the t dependence of E{Q{t)) = ( Mi (t) \Q\ Mi (t) ) 

10 




Figure 3: The same model as in figure ^ and figure < t < 1.2 s. The thin sohd hne is a graph 
of Tr{pr(t) ■ {Ir + o"3)/2), the probabihty of observing the target state | c<j ) = | (13 = +1 ) in the 
example of Grover's algorithm. The dashed line is a graph of Tr{pr{t) ■ (/^ — o"3)/2), the probability 
of observing the "undesired" output \ a^ = —1 ) . The upper and lower bounds on the probability 
of observing the target state are represented by the thick solid lines Ai(t) and \2{t)- 

and of the angle of polarization 7(t). The example of figure lU a suggests that the mean value 
of speed derived from asymptotic considerations correctly describes the average behavior of the 
"clocking" subsystem also for finite values of < t < s. As figure EJb shows, the "clocked" sub- 
system system cr is, in turn, driven, on the average, into uniform rotational motion. 




25 50 75 100 125 150 




(a) 



(b) 



Figure 4: Same parameters as in figure (a) E{Q{t)) (solid line) compared with the dashed 
straight line of slope 8/(37r). (b) The polar angle 7(t) of the Bloch vector (jTfjl as a function of t. 

Our model is so simple that we can explicitly study how the above semiclassical picture (in which 
the time parameter t acquires operational meaning from its linear relation with mean values of 
configurational observables of clocking and clocked subsystem) is distorted by a measurement per- 
formed on either subsystem. The observations made at the end of the previous section about the 
effect of reading the clock can indeed be complemented by the examination of the effect of reading 
the register. 

Suppose that the observable a^ has been measured at time r and the result +1 has been found: 
the Bloch diagram of figure Ela shows then that the evolution of the register proceeds in much 
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the same way as in the undisturbed situation of figure Q (with the only obvious difference that 
the post-measurement initial condition | 61 (r) ) lies on the unit circumference). If, instead, the 





(a) 



(b) 




1 1.2 



Figure 5: Same parameters as in figure ^ measurement of the observable a^ at time r. Frames 
(a) and (b) represent the evolution for r < t < 4r of the Bloch vector when measurement returns 
+ 1 and —1 respectively. Frame (c) represents the cumulative distribution functions of the speed 
V when the results +1 (solid thick line) and —1 (dashed line) have respectively been found; the 
solid thin line represents the c.d.f. of V in case of no measurement. 

result —1 has been found (figure Elb) the post-measurement evolution of the register is completely 

different from the unperturbed one. 

We conclude this section with an example of the insight that the time evolution of S{pr{t)) can give 

on the algorithm Ui, U2, ■ ■ ■ , Ug-i being performed by the machine. Suppose of using, instead of 

the assignment (0) of the primitive steps, Ui = U2 = ■ ■ ■ = f^s-i = ex.p{—iaa2/2), the alternative 

assignment 

(A for odd x ,^ ,. 

Ux= S n r (54) 

1^ n tor even x 
and (jH)). Figure ini gives, for this example, a full account of 



where A and B are given by 
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Kw|bi(t))|2 
1 




100 s 150 200 250 

(b) 




Si 



(c) 

Figure 6: Same parameters as in the previous figures; U^ as in ()54|) . The Bloch diagram refers to 
the time interval (0, r) needed to reach the first maximum in probabihty. 

the diffusive character [IB] of Grover's quantum search: the first maximum of the probabihty of 
finding the target state (figure IHla) is reached in correspondence of the first local minimum of 
entropy (figure IHlb): that the search has gone, before this instant, through a local maximum of 
entropy is shown with particular evidence by the Bloch diagram of figure IHlc. 

5 The role of initial conditions 

An initial condition of the form 

e 

I R{iy, ^0 ) = I R{1) ) ® J2 ^o(^)l C{x) ) (55) 

x=l 

with ipo having support in a bounded region A^ = {1, 2, . . . , e} C {1, 2, . . . , s} evolves, under (JH) as 

e-**^| i?(l); V-o ) = Y,Mt,x)\ R{x) ) ® | C{x) ), (56) 



x=l 
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where ilJo{t,x) solves, with the obvious boundary and initial conditions, the (discretized) free 
Schrodinger equation. The ensuing spreading of the wave packet leads to an increasing trend 
(with the exception of the effects of reflection at time t ^ s evidenced in figures IHlb and|2I) of the 
von Neumann entropy S{pr(t)) of the state 

s 

Pr{t) = Y,\Mt,x)f\R{x)){R{x)\ (57) 

x=l 

of the register. This is an undesirable feature because S{pr{t)) gives a lower bound on the Shannon 
entropy of the distribution of any observable of the register, for short on the uncertainty in any 
reading of the output. 

The models of the previous section where intended to show the above effect; in this section we 
devote some effort to the goal of decreasing it, by suitable choices of initial condition aimed at 
reducing the spreading of Q in the state ilJo{t,x). It is sufficient, for this purpose, to study only 
the cursor, evolving under the Hamiltonian 

A "'^ 
Ho = --^r+(x+l)r_(x) + r+(x)r_(x + 1). (58) 

x=l 

The point is to devise an initial condition ipQ which uses whatever additional finite amount A^ = 
{1, 2, . . . , e} of space resources is available as a launch pad for the cursor in an "efficient" way: this 
means both a high value of the expectation of V{iIjo) and a small value of the variance of ^(V'o) 
(we want the spreading of Q to increase at a low rate for a short time of computation). That 
both goals can be achieved is shown by examining the family of initial conditions, given by the 
eigenstates of a Hamiltonian of the form ()58|1 restricted to qubits in A^: 



x=l ^ ' 

The probability density of the speed Vk = V{ck) corresponding to each of the above states is easily 
computed from (J28|) : 

) = Ao,i)(^) • 

4 (3 - 2v'^ + cos (gf)) (sin (^))^ (sin((e + 1) arcsin(t;)))^ 



fvdv) = Vi)(^)- (60) 

.e+lJJ V Ve+1. 



TvVT^^ie + 1) {2V^ + COS (gf ) - 1) . 

The behavior of Vk is examplified by figure d We are taking there, as we will always do in this 
section for the sake of notational convenience, e to be odd 

e = 2n-l. (61) 

The examples of figure [3 clearly show the dispersive nature of the medium ()58|1 : they also show that 
increase of the mean value is accompanied by decrease of the variance (as shown by the increase 
in the steepness of the graph as k goes from 1 to n). An obvious choice for the initial state of the 
cursor emerges from the above example: 

2n-l 



2n-l /— 
x=l 



Cil) ) - I C(3) ) + I C(5) ) + . . ■ -(-l)"|C(2n-l)) ^ ^g^) 



,'n 
14 




78 .92 .97 



Figure 7: n = 5, e = 2?7, — 1 = 9. The cumulative distribution functions Fv^.{v) = Prob{Vk < v) 
corresponding to the densities (jHOl), for k going from 1 to n. The ticks on the v axes are E[Vi) < 
E{V2)<...<E{V,). 

conforming to the idea of packing the maximum number of wavelengths in the launch pad A^ = 
{l,2,...,e}, and having a, presumably easy to prepare, stationary state of the free XY chain 
localized in A^. 
The random variable Vn = V{cn) has probability density 



fv„ = I(o,i){v) 



(sin(2ra arcsin(i;)))^ 
7m{l — t>2)3/2 



(63) 



and, therefore, expectation value 

E{Vn) ^ 



t(- 

/-^ \ Ah - 



h=l 



Ah -3 Ah -I 



1-1 vf^ '- 

TT ^^ \Ah-3 Ah- 1 

h=n+l ^ 
1 



27rn 



The second moment of Vn is explicitly given by 



E{V^ 



1 . 

An 



(64) 



(65) 



The above considerations lead to the following asymptotic behavior, for large n, of the variance of 

4 -TT 



Vn. 



variyn) 



Ann 



(66) 



Equation ()64j) is a quantitative assessment of the cost in terms of space resources of achieving the 
first requisite of efficiency, namely high mean speed; similarly, (jMj) gives the cost of decreasing 
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20 40 60 80 100 120 140 

Figure 8: Solid lines: the expectation value of Q in a state | M„(t) ) evolving from an initial 
condition having the cursor in | c„ ) ; the slope of the initial linear part of the graph is correctly 
predicted by (jMj) . For comparison purposes the dashed lines show (Mi(t) \Q\ Mi(t) ) as a function 
of time and the corresponding linear fit with slope given by ()21|) . 



the variance of Vn- 

Incidentally, as the observable Q has, in the state | c„ ), expectation value 

E{Qn) = {Cn\Q\Cn) =n 



and variance 



var{Qn) 



n 



(67) 



(68) 



equation (jM|) can be read as saying that, in the initial state | c„ ), the position- velocity uncertainty 
product is given by 

var{Qr,)var{V^) ^ !l!A^ (69) 

Figures |H1 and IHl show the relevance of the above asymptotic considerations for the case of finite 
e and finite s for t < s. The effect of the initial condition is most evident if we compare the 
evolution of the state of the register from the initial state | Mi ) = | -R(l) ) | C{1) ) with the 
evolution starting from 

|M„) = |i?(l))®|c„). (70) 

This is done in figures El and El in the same probability- entropy- Block format as in figure El We 
examine there two different ways of using an additional amount s — A^ of space, of size comparable 
with the minimum amount A^ required by the algorithm. Figure ^1 summarizes the experience 
developed in |2j on the effect of using all this additional space as a telomeric chain or "landing strip" : 
as long as the cursor stays in this region the register remains acted upon by the optimal number 
of primitives. Figure ITTl shows the improvement obtained by investing part of the additional space 
as a "launch pad" on which to prepare a state in which the spreading of the cursor increases (see 
figure EI) at a lower rate than when starting from position 1. 

Comparison of figures^lc and^^c, in particular the improvement of the behavior after reflections 
at site s, shows that the idealized scenario of reversible computation (the cursor, "going back and 
forth", "does and undoes" the reversible computation) is within reach, with, as ()66|) shows, a 
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(a) 



20 40 60 80 100 120 140 

(b) 



Figure 9: s = 50, n = 5. (a) The variance of Q in a state | M„(t) ) as a function of t, compared 
with its best fit of the form const. + t^(4 — 7r)/(47rn), in the time interval (e, s — e) in which 
boundary effects can be neglected, (b) The variance of Q in the state | Mi(t) ), compared with its 
approximation ^^(3/4 - (8/(37?))^), suggested by ((22I)- 

polynomial cost in space. We note, in figure ^1 that we can do much better than in figure El 
with the same expenditure of space resources, in approximating the reversible scenario if, instead 
of the initial state (l62ll. we set the cursor in the initial state 



7n) 




2 

3n 



2n-l 



x=l 



IT 

sm I —X 
2 



C{x)). 



(71) 



The state | 7„ ) emerges quite naturally as a three-mode approximation (a linear combination of 
I c„ ) and I Cn±i )) of the initial condition that maximizes the mean speed of computation for fixed 
length e of the launch pad. 

6 Number of particles 

In the previous section we have provided examples of the benefit of spreading the initial wave 

function of the cursor {N^ = 1) on an initial launch pad instead of, as it would be classically 

"obvious", having it strictly localized at site 1. Equality (IMjl is, in this context, a quantitative 

assessment of the cost, in term of space resources, of implementing Feynman's ballistic mode of 

computation. 

In this section we abandon, in the same spirit, the classical prejudice of having a single clocking 

excitation, and present a preliminary analysis of the idea of starting the cursor in an initial state 

with A^3 > 1. The idea is to follow the motion of a swarm of several clocking agents (cursor spins 

in the "up" state) acting on the register. Stated otherwise, with reference for simplicity to the case 

A^3 = 2, we allow the clock to perform a quantum walk on the graph having the vertices (xi,X2), 

with 1 < a;i < a;2 < s, with edges between nearest neighbors [24j. 

We recall, mainly in order to establish our notation, a few elementary facts fT^ about the XY 

Hamiltonian (j58p . 

The eigenstates of Hq in the subspace N^ = n are labeled by subsets of size n of A^ = {1, 2, , . . . , s}; 

if K = {ki, k2, . . . , kn} is such a subset (where we will always assume 1 < ki < k2 < . . . < kn < s) , 
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(c) 

Figure 10: s = 50, < t < 3s, /i = 10, A^ = |_^2^/^J (the Grover-optimal number of active steps); 
Ui = U2 = . . . = Un = exp(— zacr2/2), with a and 6 chosen as in pij) and (jl^ : Ux = h for x > N; 
initial state | Mi ) = | i?(l) ) (g) | C(l) ). 



an eigenstate of Hq befonging to the eigenvalue 

n 

is given by 



i=i 



(72) 



(73) 



\E^)= Yl ViK, M)\M). 

MCAs;\M\=n 

For M = {xi, X2, . ■ ■ , Xn}, with 1 < xi < X2 < . . . < x„ < s, we have indicated above by | M ) 
the simultaneous eigenstate of T3(1), r3(2) . . . , T^i^s) in which only the spins in M are "up", and we 
have set: 

V{K,M) = det(\\v,Xx,)l,j=i,...,n) (74) 
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Figure 11: s 

Ue = f/e+l = . . . = f/e+Af-l 

|M„) = |i?(l))®|c„). 



(c) 

= 50, < t < 3s, /i = 10, A^ = 25, a and ^ as in figure CUl 
exp{-iaa2/2); U^ 



n 



5, e 



2n - 1; 



Jr, for 1 < a; < e or x > e + A^; initial state 



where the functions Vk have been defined in (fT^ . 
We set 

Wil 1^1) 3^2, • • • , 3;„| / ^ Xj| Xi, X2, • • • , Xn ) ■ I ' 5j 

It is easy to study, by the techniques of section |21 the asymptotic (as s -^ +oo and t -^ +oo) joint 
distributions of the observables Qi, and therefore to give quantitative estimates of the correlation 
between the speeds of different particles and its dependence on the initial condition. To quote just 
one example, in the subspace Ns = 2 and in the state | {1,2} ) the velocities (V^i, V2) of the two 
"up" spins (the hmits in law of Qi/t and Q2/t, respectively) have joint probability density 



fViy2iVl,V2) = I(0,V2)ivi)I(^o,v^){v2) 



Q4vlv^{2 - vf - vl) 

^v(i-^?)(i-^|y 



(76) 
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(c) 
Figure 12: Same parameters as in figure ITTl initial state | -R(l) ) ® | 7n ), as in (f7T|) . 

It is immediate from ()76p to compute the conditional expectation £'(Vi|V2) of the velocity of the 
leftmost particle given the one of the rightmost particle; it turns out to be: 



EiV^\V2) = ^ + 0{Vi). 



(77) 



In this section we advance the following idea: if the issue of the computation is the application, 
for a given number g of times, of a given primitive G to the register, initialize the cursor in the 
N^ = g subspace, in the state, say, \ {1,2, . . . ,g} ); let then the system evolve according to the 
Hamiltonian: 



A 



s-l 



H = --^U^(» r+(x + l)r_(x) + U~^ ® r+(x)r_(x + 1) 



x=l 



where 



U^ = G^^O'% for a fixed xq >g,G^ = h- 



(7J 



(79) 



An implementation of this approach is shown by the probability-entropy-Bloch diagram of fig- 
ure Uni Simple expressions for the quantities shown in figure [121 can be obtained by the explicit 
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Figure 13: fi = 4, g = 3, Xq = 6, s = 20, < t < 4s; U^, = G 
given by (jl^ ad (HI)), U^ = h for x ^ Xq] initial condition \R{1)) i 
by 



exp(— 'iao"2/2) with a and 9 
{1,2,3}) with \R{1)) given 



form of the eigenvectors of the Hamiltonian described by (|78|) and (j79j) in every eigenspace of 
A^3. For instance in the subspace A^s = 3 a complete set of eigenstates is given, for ^ = ±1 and 
1 < j < fc < s, by: 



C,]E{j,h,k}) = Yl V{{J,h,k},{xi,X2,X3})- 

1<X1<X2<X3<S 



where f? is the unit step function defined by: 

^x) = 



1, if X > 
0, if x < 



(80) 



^1) 



The spectral structure (IHUj) is peculiar of the extremely simple situation (f?^ (just one active link) 
considered there. As soon as we have more than one active link, say the primitive A acting on link 
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Figure 14: /i = 4, s = 20, < t < 4s; (a) solid line: Ua = A, Ub = B, a = 6, 6 = 8, same 
initial condition as in figure UHl (b) solid line: the N^ = 3 state has been prepared by setting the 
initial chain {1, . . . , 6} in its ground state. For comparison purpose figure IT^a is reproduced in 
both frames as a dashed line. 

(a, a + 1) and the primitive B acting on link (6, 6 + 1), with 6 > a + 1, a new phenomenon (that for 
simplicity we discuss in the N3 = 2 case) takes place: the energy eigenstates have not anymore the 
form of a linear combinations of tensors products of the form M{xi, X2) {(Js = () ®\ {xi, X2} ), with 
M{xi,X2) a monomial in A and B; related to this, the coordinates Xi, X2 lose, strictly speaking, 
the meaning of relational time ^^: given that at a given value of t, Qi = xi and Q2 = X2 we can 
only claim that the state of the register has been acted upon by a polynomial in A and B. 
This phenomenon is easily understood in terms of the Dyson expansion of the propagator: the 
probability amplitude for the two excitations being in xi, X2 (both larger than b), given that at 
time they were in yi, y2 (both < a), receives contributions not only from Feynman paths along 
which the rightmost excitation goes past a and b and then the leftmost excitation goes past a 
and b (along such a computational path the state of the register is modified by BAB A), but also, 
among others, from paths along which both excitations go past a before both going past b (along 
such a computational path the state of the register is modified by BBAA). 

Waiting for an algorithm that might benefit from the above possibility of simultaneously exploring 
different computational paths (concurrency ?), we explore, in figure El the idea (or classical 
prejudice?) that this nuisance can be in part avoided by using suitable initial conditions. The 
idea, suggested by (jUUIl . is of course to prepare on A^ = {1, 2, . . . , a} an initial N3 = g state such 
that the excitations travel as spatially well localized wave packets of so different speeds that it is 
at any time unlikely that they simultaneously hit the region (a + 1, 6). 

We conclude with a remark about our insistence, throughout the paper, in gathering experience 
about the behavior of the evolution of a state of an initial subchain A^ = {1, 2, . . . , e}. 
We observe that an initial state (not necessarily in the N^ = 1 subspace) in A^ of the form 



in ) 




z£{-i,iY jeM 



T,{X) = (-l)^-(-) ), 



(82) 



where Im is the indicator function of the set M and x = 1, . . . ,e, can be prepared as a post- 
kickback state (with respect to an ancilla qubit) after the reversible evaluation of a function 
/ : {—1,1}'' — > { — 1,1}. We conjecture that subsequent evolution oi \in) under the Hamiltonian 
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(IHF^ on As = {1, 2, . . . , s}, with s >> e, might help in setting tests of hypotheses about the Fourier 
coefficients 

z6{-l,l}^ jeM 

of the function / via time-of-flight techniques. There is at least one non trivial case in which 
the above conjecture works: having prepared all spins in {e + 1, . . . , s} in the "up" state, the 
Deutsch-Josza alternative ^9] "constant (c0 = 1) vs. balanced (cg = 0)" becomes equivalent to the 
alternative "stationary vs. non stationary" under the Hamiltonian (|3Hj). about the state of the 
overall system. 

7 Conclusions and outlook 

The pure XY Hamiltonian Hq given in ()58|) describes, in the Luther-Liischer-Susskind formalism 
)221 1211 1^ . a massless Dirac quantum field on a 1-dimensional lattice. The full Hamiltonian (0) 
is suggestive of the minimal coupling of this Fermi field, implementing the clock, with additional 
quantum fields implementing the register. This work intends to contribute to the line of research, 
that seems to be emerging these days [HI I2HI 1^0] , devoted to making this connection between quan- 
tum computing and relativistic quantum field theory explicit. It is an easy guess that this quantum 
field theoretical intuition was well present in the original work ^Uj- Particularly penetrating is, 
in this respect, Peres' remark that in Feynman's model calculations run forward and backward in 
time just as particles and antiparticles in Feynman's classical work on relativistic quantum field 
theory {2Q, p. 3269). As a further remark, we observe that the three-body interactions needed 
by Feynman's model are hard to conceive out of a field theoretical context. 

It is because of this field theoretical perspective that we have tried to avoid any "engineering" (space 
dependence) of the coupling constant A in (jHHl), well aware of the fact that, in the Dirac — > XY 
correspondence, A is related to the spacing adopted in the lattice approximation. In such a context 
it would be very hard to understand (without a projection mechanism [7], which seems to have an 
exponential cost) the implementation of a space dependence such as 

A(x) = const. \/x{x — s) (84) 

that leads in j^Hl and [7j to the existence of sharply distinguished instants in which the position 
of the cursor is certain. Nor would it be easy to understand ()84|) in a solid state implementation 
j3], where A is related to the effective mass of the clocking excitation. 

In this paper we have focused our attention on the clocking field t{x), singled out as the one 
which, under suitable boundary conditions and for initial conditions localized close to the bound- 
ary, exhibits particle-like excitations performing, for long enough intervals of t, a quantum walk in 
a distinguished direction. 

Spatial homogeneity of the chain leads to the existence of the limit in law V = limt^+oo <5(t)/t 
for the position of such an excitation on a semi-infinite {s — > +cxd) box. In the N^ = 1 subspace, 
because of Peres' conservation law j^ni; the observable Q acquires the meaning of relational time 
(given the observed value of Q, the state of the register is known with certainty) and, therefore, 
the random variable V acquires the meaning of number of computational steps per unit t. The fact 
that the variance of V is strictly positive has the effect that in terms of the parameter time t (as 
opposed to relational time Q) the evolution of the register appears to be dissipative: we have, for 
a simple model, written the corresponding Lindblad evolution and studied the ensuing build-up of 
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entropy. 

On our simple instance of quantum search we have shown that, in the "low level", physical ap- 
proach that we pursue (in which time runs, for the cursor, because it is coupled with an additional 
quantum field) the build-up of entropy imposes an upper bound on the probability of finding the 
target state which is more severe than the one predicted by the "high level" , algorithmic approach 
(in which the successive primitives are applied by an external macroscopic agent). 
In the attempt of decreasing the deficit in the probability of success in a quantum search, due to 
the decohering effect of the coupling with the clocking field, we have provided examples of the 
benefit of spreading the initial wave function of the cursor on an initial launch pad instead of, as 
a classical prejudice would suggest, having it strictly localized at one site. 

We have, similarly, abandoned the classical prejudice of having a single clocking excitation, pro- 
viding a preliminary analysis of the idea of starting the cursor in an initial state with A^3 > 1. We 
have shown, in this context, an efficient way of iterating the application of a single primitive to the 
register and experienced the possibility, by a suitable choice of the initial conditions, of steering 
the quantum walk of the excitations in such a way as to reduce the conflicts about the order of 
application of non commuting primitives. 

The case N3 > 1 deserves, we think, further research, both from the algorithmic and the physical 
point of view. 

From the algorithmic point of view we plan to examine other instances (beyond the one cursorily 
examined at the end of section 6) in which time-of-fiight spectroscopy (based on the Fourier trans- 
form vs. speed relationship recalled in section 3) of the post-kickback state can answer Yes/No 
questions about the algorithm. 

From the physical point of view, the "obvious" choice of the "all down" reference state made 
throughout the paper is far from being optimal from the point of view of studying the thermody- 
namic cost of resetting the register. The best reference state for the study of this ultimate cost of 
reversible computation would of course be the ground state and, for Hamiltonians of the form (j3Hj) . 
with s even, it is an N^ = s/2 state. This will require, we think, the formulation of an appropriate 
Bethe Anzatz for the Hamiltonian (jnj. 
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